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The effects of an imposed, axial magnetic field Bzo on hydrodynamics and energetic 
electrons in inertial confinement fusion (IGF) indirect-drive hohlraums are studied. We 
present simulations from the radiation-hydrodynamics code HYDRA of a low-adiabat 
ignition design for the National Ignition Facility (NIF), with and without Bzo = 70 
Tesla. The field’s main hydrodynamic effect is to significantly reduce electron thermal 
conduction perpendicular to the field. This results in hotter and less dense plasma on 
the equator between the capsule and hohhaum wall. The inner laser beams experience 
less inverse bremsstrahlung absorption before reaching the wall. The x-ray drive is thus 
stronger from the equator with the imposed field. We study superthermal, or “hot,” 
electron dynamics with the particle-in-cell code ZUMA, using plasma conditions from 
HYDRA. During the early -time laser picket, hot electrons based on two-plasmon decay 
in the laser entrance hole ( Rngan et al\\201(h are guided to the capsule by a 70 T field. 
12x more energy deposits in the deuterium-tritium (DT) fuel. For plasma conditions 
early in peak laser power, we present mono-energetic test-case studies with ZUMA as 
well as sources based on inner-beam stimulated Raman scattering. The effect of the field 
on DT deposition depends strongly on the source location, namely whether hot electrons 
are generated on field lines that connect to the capsule. 


1. Introduction 

Using a magnetic field t o enhance inertial fu sion is an old idea ( Jones fc MeadI 1986h 
receiving renewed interest ( Slutz fc Vesevl[^Q12[ ). An imposed field is being i nvestigated at 


LLNL as a way to imp rove capsule performance and achieve ignition on NIF (|Perkins et al 


20131 . 2Q14 Ha 2Q15 L These simulation studies show an initial field of 40-70 T increases 


both the likelihood of ignition and the fusion yield by reducing electron-heat and alpha- 
particle loss from the hot spot. Earlier experiments at the Omega laser facility with 
an imposed 8 T axial held show increased fusion yield and ion temperature in spher¬ 
ical implosions ( Chang et a/.ll201ll: iHohenberger et a/]l2012h . The held may also limit 
hydrodynamic (e.g. Rayleigh-Taylor) instability growth, and reduce the negative effects 
of the growth that does occur. The held also increases the plasma temperature in the 
unde rdense hohhaum hh. whic h could reduce SRS and improve laser propagation to the 
wall (|Montgomerv et a/.N2Q15h . A pnlsed-power appr oach is being developed to impose 

and is sketched in Fig. [H Laser- 


Bzo = 70 T on a NIF hohhaum 
driven capacitor-coil syste ms are a possible way to impose 100-1000 T helds (|Fnjioka et al 
2ni.^lPoi]ock et fl/Jlimlh . 
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Figure 1. Sketch of pulsed-power coil design to impose axial magnetic field on NIF hohlraum. 
Left: hohlraum surrounded by solenoidal coil, helded on a Diagnostic Instrument Manipulator 
(DIM). Additional hardware needed for holding has been removed. Right: diagram of hohlraum 
and solenoid, with full holding hardware included. The red region indicates the gold or other 
high-Z hohlraum wall, while the gray regions outside the coil are additional support structure. 


This paper presents simulation studies of how an imposed held affects hohlraum hy¬ 
drodynamics and energetic elec trons. First, we repor t on simulations using the radiation- 
hydrodynamics code HYDRA ( Marinak et al 200 ih with and without an imposed held 
of igniti on experiment N120321. Then we sh ow studies with the particle-in-cell code 
ZUMA ( Larson et al 2Q1QI: Strozzi et al 2012 ) of the held’s effect on energetic or “hot” 
electrons. 

We study NIF shot N120321, which used a 4-shock, low-adiabat or “low-foot” laser 
pulse, a plastic ablator, and a cryogenic DT ice layer. It achieved the highest fuel areal 
densit y to date on NIF, and has been extensively modeled to understand its low neutron 
yield dciark et a/.1l2Q15 ). Here, we use hydra’s magnetohydrodynamics (MHD) package 
( Koning et a/.ll2006l) with = 70 T, which was not present in the actual experiment. 
We include the J x B magnetic pressure force, a simple Ohm’s law E = r]J — v x B^ 
Ohmic heating, and anisotropic electron thermal conductivities parallel and perpendicu¬ 
lar to B (but not the Righi-Leduc heat flow along B xVTq). This neglects several effects 
which could be important, and will be studied in future work, namely the self-generated 
or “Biermann battery” dtB oc VTe x Vrie field and the Nernst effect {E (x B x VTe). In 
our runs, the B field roughly follows the MHD “frozen-in law” for the highly-conducting 
plasma flow. The primary effect of the field is to reduce electron heat conduction per¬ 
pendicular to B. This leads to a hotter hohlraum fill, and a wider channel between the 
capsule and hohlraum equator. The inner cone of beams (pointed toward the equator) 
better propagate to the wall, which gives more equatorial x-ray drive and a less oblate 
imploded capsule. This would reduce the need for energy transfer to the inners, and 
probably reduce their backscatter - both due to the lower power and higher temperature. 

Besides hydrodynamics, we also study hot electron dynamics. Hot electrons are a 
generic aspect of intense laser-plasma interactions (LPI). They are produced in any para¬ 
metric process that drives a Langmuir wave. Of particular interest in IGF are stimulated 
Raman scattering (SRS) and two-plasmon decay (TPD). These are the decay of a light 
wave to a Langmuir wave and, respectively, a scattered light wave (SRS) or another 
Langmuir wave (TPD). In many laser-produced plasmas, the daughter Langmuir waves 
are damped primarily by collisionless Landau damping, which entails the resonant inter¬ 
action of the wave with electrons at its phase velocity. This is typically greater than the 
electron thermal speed, and therefore produces a population of superthermal or “hot” 
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electrons. Experiments show the resulting hot-electron spectrum from a single paramet¬ 
ric process is roughly exponential with “temperature” dN/dE oc g{E)e~^l'^^ {E is 
the hot electron kinetic energy), with g = for a non-relativistic Maxwellian. NIF 

experiments with gas-filled hohlraums have shown hard x-ray output consistent with a 
two-temperature hot-electron pop ulation, attributed t o Raman backscatter, and TPD 
or SRS at quarter-critical density (|DQDPner et a/.ll2Q12h . Relativistic processes that pro¬ 


duce > MeV electrons at intensities IX^ > 10^^ W cm“^ are of great interest in the 
short-pulse and fast-ignition fields, but are not discussed here. 

This paper focuses on hot electrons in ignition hohlraums, though similar considera¬ 
tions apply to directly-driven targets. Hot electrons impede ICE in several ways, namely 
implosion asymmetry and fuel preheat. The laser power transferred to hot electrons 
generally stays in the target, so is not a power loss like backscattered light. But, the de¬ 
position in space and time differs from the intended inverse-bremsstrahlung absorption 
of the incident laser. NIF hohlraums with high hohlraum gas fill density (> 0.9 mgcm“^ 
He) have generally shown large SRS from the inner beams. This reduces the inner-beam 
power reaching the wall - both by scattered light and Langmuir waves - which makes the 
implosion more oblate (or “pancaked”). The Langmuir wave energy remains in the target, 
but heats the hohlraum wall by conduction in a much larger area than the inner-beam 
spots. To control symmetry, cross-beam energy transfer has been used to move power 
from the outer to inner beams, inside the target (|Michel et a/.ll2QQ9l ). Hot electrons with 
energy > 170 keV can also preheat the fusion fuel (e.g., cryogenic DT ice layer) by 
depositing energy separate from the inten ded shock sequence and capsule compression 
(|Salmonson et a/.ll2010l: iHaan 12011 ). This results in a higher fuel adiabat, which 
significantly reduces the achievable compression. 

We propagate hot electrons with the hybrid-PIC code ZUMA through plasma condi¬ 
tions from HYDRA. We run ZUMA in a “Monte-Carlo mode” with no E or H fields, 
except sometimes a static B. Hot electrons undergo energy loss and angular scattering 
as they propagate, and the energy deposition profile is found with and without an initial 
Bzq. We first present an unphysical test-case study of mono-energetic hot electrons di¬ 
rectly incident on the capsule (unrealistic for LPI-produced hot electrons) early in peak 
laser power (time 18 ns). A minimum initial energy Eq =125 keV is needed to penetrate 
the ablator and reach the DT layer. The maximum energy deposited in the DT layer, 
Edt^ occurs for Eq =185 keV and is Edt/Eq = 13%. Higher energy electrons do not 
fully stop in DT. 

We then examine a realistic hot-electron source, consistent with two-plasmon decay 
during the early-time “picket” or initial part of the laser pulse (time 1 ns). The deposition 
is mostly in the high-Z wall, as expected from solid-angle arguments, and Eor/Eh = 
2.2 X 10“^ with Eh the total injected hot electron energy. Adding a uniform 70 T axial B 
field strongly magnetizes the hot electrons in the hohlraum fill, guides them to the capsule, 
and increases Eor/Eh by 12x to 0.026. This may not degrade fusion performance, since 
NIF experiments have shown greatly reduced picket hot electrons with pulse shaping, e .g. 
a low-power “toe” to burn down the window ( Dewald et ~^l2015l : iMoodv et a/.l 12014 ). 

Finally, we consider a hot-electron source consistent with SRS of the inner laser beams, 
early in peak laser power (18 ns). With no the hot electrons deposit throughout the 
target, with a very small EjjT/Eh ~ 1.2 x 10“^. With Bzo = 70 T, the field strongly 
magnetizes the hot electrons in the hohlraum fill gas. The deposition in DT is greatly 
increased (decreased) for hot electrons originating on field lines that do (do not) connect 
to the capsule at this time. 

This paper is organized as follows. We describe our MHD simulation methodology in 
section [2] and our MHD results in section [H The ZUMA simulation method is detailed in 

















4 


D. J. Strozzi 




Figure 2. Incident laser power on NIF shot N120321 on the inner (dashed) and outer (solid) 
laser beams during the early-time picket (left) and peak power (right). The dotted vertical 
lines at 1 and 18 ns indicate times HYDRA plasma conditions are used for ZUMA hot electron 
studies. 


section HI Section [5] discusses test cases of mono-energetic electron propagation through 
the capsule. Section [ 6 ]presents ZUMA results for a TPD-relevant source during the picket, 
and shows a 12x increase in Ejjt with a 70 T axial field. In section [3 we present ZUMA 
results early in peak laser power with an SRS-relevant hot electron source, using plasma 
conditions and the B field from our HYDRA simulations, and find a strong dependence 
in Edt on source location. We conclude in section [HI 


2. HYDRA MHD simulation method 

We use the radiation-hydrodynamics code HYDRA to simulate the NIF hohlraum 
experiment N120321. This shot used a 4-shock, “low-foot” laser pulse (shown in Fig. 
El, a plastic ablator (C 0 . 42 H 0.57 plus small amounts of O impurity, and Si dopant to 
control x-ray preheat) with a DT ice layer, and a depleted uranium (DU) hohlraum with 
a thin 0.7 /im inner gold coating. The hohlraum fill gas was 0.96 mgcm“^of He. The 


meth odology is the standard one in use for hohlraum simulations at LLNL (| Jones et all 
I 2 ni 2 h . entailing the “high-flux model” with detailed configuration accounting (DCA) 
for non-LTE ma terial properties, and an electron thermal flux limit of 0.15 the free- 
streaming value ( Rosen et al ] l201ll) . The runs use a 3D mesh with one zone and periodic 
boundary conditions in azimuth, and are effectively cylindrical 2D (r-z). The mesh is 
managed with an arbitrary Lagrangian-Eulerian (ALE) approach, designed to keep the 
mesh as Lagrangian as possible. We use the full, incident laser energy of 1.52 MJ, and 
neither remove measured backscatter n or degrade the las er power as is needed to match 
x-ray drive data in high gas fill targets (| Jones et a/.ll2012f ). Our HYDRA simulations are 
therefore not proper post-shots, but address the role of an imposed field in ignition-scale 
designs, and provide relevant plasma conditions for hot electron studies. 

A distinct aspect of the present work is the inclusion of an initial axial mag netic field 
Bzo. hydra’s magnetohydrodynamics (MHD) package (|Koning et was used 

to model Bz and the resulting Br required by V • H = 0. No azimuthal field is produced 
for our axisymmetric geometry and simple Ohm’s law. The MHD package uses a 3D finite 
element method, with appropriate boundary conditions to be effectively axisymmetric. 
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We used the Ohm’s law 

E = r]J — vxB (2.1) 

where 77 , J, and v are scalar resistivity, net current, and center-of-mass velocity, respec¬ 
tively. The B field is evolved via dtB = —V x E. The MHD package as used here affects 
the matter in three ways: a) the J x B force, b) r]J^ (Ohmic) and other heating terms, 
and c) a tensor electron thermal conductivity: 

/ ^ B 

K, = y — bbj K,\\bb b = (2.2) 

The Righi-Leduc effect, with a separate conductivity along BxVTq^ is currently neglected(Te 
is electron temperature). An artificial flux limit / is imposed, as is typical in hohlraum 
simulations. Specifically, the component of the heat flux pe along each logical index (not 
physical) coordinate i is limited: qe'i = min[i • k • VTe, fqrs] where qps = neTe^‘^ 
is the free-streaming heat flux. 

The anisotropic heat conduction has the largest effect in our simulations. We expect 
5 > 1 T to significantly reduce below its unmagnetized value The Hall parameter 
H = uJceHi for thermal electrons is 


H — uJqqTqi — 


B 


Bo = 


^eff = 


(327r)^/^ He 

3 (47reo)^ Te 

EJ,Zf\nAei 


(2.3a) 

(2.36) 


(2.3c) 


For each ion species Zi is the ionic (not nuclear) charge, Ui = /^n/, and nj = 
is the total ion number density. In practical units, Ho[T] = 4.73(ne/ncr)^eff/Fe[keV]^/^ 
with Ucr = 9.05 X 10^^ cm“^ the critical density for light of wavelength 351 nm. For 
Zi = 2 Be at Uq = O.lncr and Tg = 3 keV, typical of the underdense hohlraum fill, we 
find InAei = 7.9 and Bq = 1.43 T. Given Bzo = 70 T, most of the underdense plasma 
fill should be strongly magnetized. decreases with H according to 


K,± ^ 1 + PiH 

1 + P 2 H + P3H‘^ + p^H^ ’ 


(2.4) 


with the Z^-dependent fitting coefficients pj given in lEpperlein fc Hain"^ (Il986l). For 
Zi = 2 He, = 0.1 for H = 1.6. n\\ is found either from the Lee fc More |l984l) 

formulation to include dense-plasma effects, or interpolation from an advanced table. 
The Epperlein and Haines results are used to include electron self-collisions {Zi < 00 ) 
and dependence on H. 


3. MHD simulation results 

The HYDRA runs of NIE shot N120321 with and without MHD are qualitatively 
similar. The principal difference is the MHD run has higher electron temperature in 
some regions and a wider channel of He fill gas at the equator. This leads to better inner 
beam propagation to the wall (less inverse-bremsstrahlung absorption in the low-Z fill), 
and results in a less oblate capsule. Eigure [3] shows the material regions and electron 
density for the two runs at 18 ns, during the rise to peak power. We use this time 
for ZUMA simulations of SRS hot electrons. A more detailed density plot is in Eig. 0 
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Figure 3. Plasma conditions from HYDRA simulations of NIF shot N120321 at 18 ns, used 
in SRSPEAK ZUMA run series. The left {z < 0) half of the hohlraum is plotted, but the 
simulation included both halves. Top half (r > 0): without MHD, bottom half (r < 0): with 
MHD and initial axial magnetic field Bzo = 70 T. Left: material region: DT, He, CH (two green 
regions). An, DU label deuterium-tritium, helium, plastic, gold, and depleted uranium. Right: 
free electron density in units of critical density for 351 nm light. 


Electron temperature Te with and without MHD at several times during peak power is 
displayed in Fig. El The field increases Te mainly in the laser-heated gold, such as the 
outer-beam “bubble” (r, z) « (0.2, —0.3) cm. The He gas fill is also hotter with the field, 
but less so than the gold. The low density plasma outside the laser entrance hole (LEH) 
is cooler with the field, though the total energy in this region is small. 

The B field for the MHD run is plotted in Fig. [5l It roughly follows the MHD frozen-in 
law, and advects with the radial motion of the ablator and high-Z wall. The compressed 
field approaches 300 T and continues to grow with time. The white stream line that just 
touches the capsule outer radius at (z^r) = (0,0.1) roughly separates field lines that are 
still connected to the capsule (r < 0.1 cm at ^ = 0), from those that have advected with 
the ablated blowoff and no longer connect to the capsule. 

The capsule density at the end of peak power is plotted in Fig. [6l The dense fuel is 
oblate without the field, but becomes close to round with it. 


4. ZUMA Hot electron simulation method 

We propagate hot electrons through fixed plasma conditions from HYDRA using the 
hybrid-PIC code ZUMA in a “Monte-Carlo” mode. We do not include forces from E 
and B fields, except when we include a specified (static) B field. The background plasma 
properties are not updated. In other work, HYDR A and ZUMA have been coupled to run 
in tandem, and applied to fast ignition designs (|Strozzi et a/.l[2Ql^ . The hot electrons 
undergo collisional energy loss off background electrons, and angular scattering off back¬ 
ground electrons and ions. We neglect collisions among hot electr ons, since their densit; 


is much less than the background species. We use the formulas in lRubinson et al\ (|2Q1 
for a fast electron with v vtc- The energy loss rate (stopping power) is given by 


% 


dE CqTIq 

Eu ^ -Td 

dt rrieV 


(4.1a) 
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Figure 4. Left: electron temperature in keV from HYDRA MHD simulations shown in Fig.fS] at 
several times. Right: temperature difference with MHD minus without MHD. Light blue contour 
marks region boundary for DT, CH, and An. Dashed contours are temperature differences of 
±0.5 keV. 
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Figure 5. \B\ at 18 ns from HYDRA MHD simulation shown in Fig. [S] The dark-red color in 
the capsule is the uncompressed held of 70 T. White curves are stream lines (integral curves) 
of the vector held {Bz,Br). Light blue boxes and text indicate ZUMA hot electron sources for 
runs with and without MHD. 



Figure 6. Capsule shell density at 21.5 ns from HYDRA simulations shown in Fig. [3] without 
MHD (top half, r > 0), and with MHD and Bzo = 70 T (bottom half, r < 0). The shell is 
oblate without MHD, while it is close to round with MHD. This rehects the improved inner 
beam propagation with the held. 


CqTIq 


:ln- 


E 


Ld = In 


yj2meE hzOp ’ 

pv In 2 

hujp + 1 2 


hujp E ^ rrieC^ 

9 ln2 + 1/8 f 1 


+ 16 + 


- 1 

27 


(4.1 &) 
(4.1c) 


Ce = e^/dTrcg, the fast electron kinetic energy E = mgC^e, 7 = e +1 is the Lorentz factor, 
and p = 'yrUeV. rie is the total (free plus bound) background electron density. Ld given 
above is valid for energy loss off free electrons, or bound electrons for sufhciently high 
E or Ue (the “density ehFect”). This assumption may not be valid for all electrons. The 
angular scattering rate is 


d{e^) 


2Ce 

p‘^v 


^Tii ^ Lqi T tIq^^L 


(4.2 a) 


dt 
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Lsi = in — -— u.Doy'i;“/c“, 



< E < meC^,(4.26) 
(4.2c) 



(4.2 d) 


Uej is the free electron density, (Z^) = fiZf and we use the same ion species notation 
as after Eg. 12.31 except Zi is the nuclear (not ionic) charge. Ig is a screening length, which 
we take to be the free electron Debye length. For neutral atoms, it should be replaced 
by the atomic radius. In any event, we impose a minimum of 1 on Lgi^ and Lge- 
We run ZUMA in 2D cylindrical geometry. ZUMA currently operates with constant 
(but different) grid spacings dr and dz. The HYDRA plasma conditions a re interpolate d 
onto a uniform mesh with dr = dz = ^ /am using the OVERLINK package ( Grandvl[l999l) . 
This small spacing is needed to resolve small features, such as the gold wall and DT layer. 
The ZUMA time step is 1 fs, which is chosen to adequately resolve the dependence of 
dE/dt on E for small E and high rie- ZUMA stops following electrons when E < 5.11 
keV and locally deposits their kinetic energy. In this paper, ZUMA injects hot electrons 
from a distribution that is a product of an energy spectrum dN/dE times a polar angle 
spectrum dN/dCl. For a thermal spectrum with a “temperature” T/j,, we use a relativistic 
Maxwell-Jiittner distribution: 


He 


C[l + e/2]^/^ [1 + e] 


(4.3) 


C is a normalization constant, and the two bracketed factors are absent for a non- 
relativistic Maxwellian. 

5. Mono-energetic electron propagation through capsule at peak 
power: CAPTEST series 

This section considers the propagation of electrons directly incident on the capsule 
during peak laser power, as a function of electron energy. We call this the CAPTEST 
series of ZUMA runs, and stress this source is not realistic for LPI-generated hot electrons. 
Rather, our purpose is to understand where electrons that reach the capsule deposit their 
energy, and which energies pose the greatest preheat risk. We use plasma conditions from 
our HYDRA simulation of NIF shot N120321 with no MHD at time 18 ns. The same 
conditions are used in the SRS-relevant SRSPEAK series discussed below in section [71 
The time 18 ns is during the rise to peak power (see Fig. and has significant inner- 
beam SRS. An analogous time in shock-timing (“k eyhole”) shots has been identified as 
possibly having a large hot-electron preheat effect (|R.obev et a/.ll2Ql^ . 

Figure [7| shows the total (free plus atomically bound) electron density in the HYDRA 
simulation. Mono-energetic hot electrons are injected in a cylinder of radius 500 /am at 
2 : = —0.1 cm, with an initial velocity in the 2 : direction. The hot electrons experience 
energy loss and (in some runs) angular scatter, but no forces from E or B fields. The 
resulting energy deposited per volume, zoomed on the capsule, is plotted in Fig. |8l The 
case with angular scattering shows large spreading of the hot electrons in the dense CH 
ablator. Since the absolute number of hot electrons introduced is arbitrary, we express 
the deposition as energy density per injected hot electron energy. 

We plot the fraction of injected energy that is deposited in various regions, or escapes 
to the boundaries, in Fig. |9l The deposition in DT (essentially the dense fuel layer, not 
the less-dense proto-hotspot) is shown in Fig. [TOl All electrons with F < 125 keV stop 
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rig total [cm 18 ns 



Figure 7. Electron density (free plus bound) at 18 ns in HYDRA simulation of NIF shot 
N120321 with no MHD (discussed in Sec. [3]), used for CAPTEST series of ZUMA runs. Bottom 
panel zooms on capsule. Material regions are separated by solid lines and labeled with text as 
in Fig. O The ZUMA hot electron source is indicated. 


in the CH ablator. This sets a minimum energy hot electrons must have when they reaeh 
the eapsule (at the time 18 ns) to reach the DT layer. Above this energy, the fraction 
deposited in DT increases, until the hot electrons have enough energy to not stop in 
the DT. With no angular scatter (left panel in Fig. [9]), hot electrons eventually cross 
the capsule, and exit the problem through the LEH. Angular scatter lowers and spreads 
out the peak in coupling to DT. It also causes some hot electrons to reach the Au/U 
hohlraum wall. In both cases, the deposition in the He hohlraum fill gas is negligible. 

A simple model illustrates the basic features of Fig.[T0l especially for no angular scatter. 
Imagine a hot electron starting at position zq in the CH ablator, with initial energy Eq 
and Vz > 0. We use a ID slab geometry with CH from zq to zi, DT from zi to 2 ^ 2 , and 
CH for z > Z 2 . We seek the fraction of initial energy deposited in DT, (pjjT = EdtIEq^ 
where Edt is the energy deposited in DT. The hot electron loses energy as it moves to 
increasing z according to dEjdz = —f/2E, where / is a constant, and we include only 
the leading-order dependence of stopping power on energy, for E rUeC^. Integrating 
from Za to Zh gives E^ = E^ — f{zb — Za). An electron with Eq < Eqi fully stops in the 
CH with z < zi, one with Eq > E 02 crosses the DT, i.e. stops at z > ^ 2 , and one with 
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Figure 8. Energy deposited by 175 keV mono-energetic hot electron source, without (top) 
and with (bottom) angular scattering, for CAPTEST ZUMA run series. 
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Eigure 9. Eraction of injected hot electron energy deposited in different regions, or escaped 
to boundaries, for CAPTEST ZUMA run series. Plot on (left, right) is (without, with) angular 
scattering. 


Eqi < Eq < Eq 2 stops in the DT layer, i.e. zi < z < Z 2 . A straightforward calculation 
gives 

0 , Eq < Eqi 

ct>DT = < [1 - Eoi <Eo< Eo2 (5.1) 

^ [1 - - [1 - Eo 2 < Eo- 

Eqi = fcH{zi - ^o) and Eq 2 = r^oi + fDT{z 2 - zi). For Eq = F^oi + 5E with 5E small, 
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Figure 10. Solid curves: (I)dt — fraction of injected hot electron energy deposited in DT, with 
(black) and without (red) angular scatter, for CAPTEST ZUMA run series. Closed circles are 
ZUMA simulation points. Dashed curves: power-law fits to E ^ 200 keV points. Thin blue curve 
is approximate form explained in text. 


4>dt ~ {25E/ , and for Eq > Eq 2 , 4>dt ~ /dt(^2 — ^i)/2F^o- The simple model for 
4>dt^ with Eqi = 135 keV and F^o 2 = 155 keV, is plotted as the solid blue curve in Fig.fTOl 
The model is close to the red, no-scattering result, though the capsule curvature smears 
the peak compared to the simple model. Figure [TO] also includes least-square power-law 
fits to the Eq ^ 200 keV results: (j)DT = (^o/52.3 keV)“^‘^^ without angular scatter, and 
0DT = (^o/h2.7 keV)“^‘^^ with scatter. These are both close to the Eq‘^ scaling of our 
simple model for Eq ^ F^o 2 - 

We apply our mono-energetic results to a thermal spectrum in Figs. [TT] and [121 
find DT preheat comes mainly from hot electrons with energies ^160 keV, for > 20 
keV. Figure Hi] shows the coupling to DT of a thermal. Maxwell-Jiittner spectrum with 
Th = 50 keV. The black curve is the DT coupling fraction from Fig. [10] (black curve there 
too), and the blue curve is the thermal energy spectrum E * dN/dE for = 50 keV. 
The red curve is their product, namely the energy coupled to DT by electrons of a given 
energy, in a thermal spectrum. The red curve exhibits behavior akin to the “Gamow 
peak” in fusion reactions, with a location determined essentially by the steeper black 
curve. Figure [12] shows the overall c^dt integrated over the thermal spectrum vs. Th. 
This peaks slightly above 5% near Th = 90 keV. The red curve is, as a function of T/^, 
the hot electron energy of maximum Edt^ i-e. the energy of the peak in the red curve 
in Fig. [TTJ This increases slowly with T/j,, and is at > 160 keV for all Th of interest. It is 
thus important to correctly model these hot electrons to calculate DT preheat, even for 
Th < 160 keV. 


6. Hot electron propagation in early-time picket: PICKET series 

This section studies hot electron dynamics during the initial laser “picket,” and the 
effects of an axial B field. The principal way hot electrons are produced during the 
picket is LPI in the LEH. This can be two-p lasmon decay (TPD ) for Ue ~ Raman 

scattering, or a multi-beam variant of it ( Michel et all [ 2 QI 5 I ) . NIF experiments have 
shown the picket hot electrons can be reduced by shaping the picket pulse, for instance 
by turning the inner beams on before the outers to blown down the window at low 
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Figure 11. Coupling to DT of a Maxwell-Jiittner hot electron source. Blue: E * dN/dE for 
Jiittner source with Th = 50 keV. Black: (^dt = energy fraction deposited in DT, with angular 
scatter (solid black curve from Fig. IIQD . Red: product of black and blue curves. 
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Figure 12. Black (left y-axis): injected hot electron energy fraction coupled to DT, for 
Maxwell-Jiittner source with temperature That- Red (right y-axis): hot electron energy with 
peak coupling to DT (i.e., peak of red curve from Fig. E]). 


power. Exp eriments at the Om ega laser studied hot electrons from TPD during window 
burn-down ( Rugan et 

The DT fuel is particularly sensitive to hot electrons produced during the early time 
picket pulse: Aentropy =Ei:)T/temperature, so a small Edt added when the fuel is cold 
produces a large entropy increase. In addition, melting the cryogenic DT layer before the 
first shock arrives causes t he inside surface to expand, which can degrade the ability to 
shock-time ( Thomasll2Q15 ). For indirect-drive ignition designs, this occurs for Edt ^0.1 
J. NIF ignition-relevant hohlraum experiments show total hot-electron energies E^ ^ 1 J 
with Th ^ 80 keV. Calculations typically show Edt/E h ^ (2 — 5) x 10“^, giving preheat 
Edt ^ (2 — 5) X 10“^ J well below melt. 

We use HYDRA plasma conditions at 1 ns, shortly after the outer-beam power has 
peaked, for ZUMA calculations. We call this the PICKET run series. Eigure [13] shows the 
material regions and laser intensity. We source the hot electrons in a 500 fim radius circle 
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Figure 13. Material region (top) and laser intensity summed over all beams (bottom) at 1 ns 
from HYDRA simulation with no MHD of NIF shot N120321, used for PICKET series of ZUMA 
runs. Hot electron source is indicated as white box, and is placed only in the z < 0 half of the 
two-sided {z < 0 and z > 0) ZUMA domain. 


at the left-side LEH {z = —0.45 cm), which is roughly the extent of high laser intensity. 
Since TPD does not generally produce collimated hot electrons, we use an isotropic 
source with velocity-space dN/dft constant {ft is solid angle in velocity) for polar angles 
between 0 and 90°, and zero otherwise (i.e., uniform in the forward-going half-space). 
The energy spectrum is a Maxwell-Jiittner with = 80 keV, which is consistent with 
hard x-ray data on NIF (discussed below). 

We expect a 70 T axial field to strongly magnetize the hot electrons in the low-density 
hohlraum gas fill, and guide them to the capsule. Recall that we inject a divergent hot 
electron source, so the question is whether the field confines them in space. It will not 
collimate them, i.e. reduce their velocity-space divergence. The electron Larmor radius 
^Le = which for E = 100 keV and 5 = 70 T is r^e = 16.0 jum-sma {a is the 

angle between B and p). This is much less than the relevant plasma scale lengths. Also, 
the cyclotron period Tee = ^ir^meleB is 0.510 ps-y for 5 = 70 T, which is much shorter 
than the propagation time through the hohlraum. Figure [TU plots Tee and the time for 
90° root-mean-square angular scatter, Tas- {0‘^) = (90°)^ for dt = Tas in Eq. 14.21 We 
consider two fully-ionized cases: one representative of the hohlraum fill: 0.96 mgcm“^ 
of He at Te=l keV, and one of the ablator: 1 gcm“^ of CiHi at 200 eV. The plot shows 
all hot electrons are magnetized in the He, while those with E > 300 keV in the CH are. 
Even if Tee exceeds Tas, r^e is much smaller than typical capsule dimensions ^ lOO’s pm. 

The energy deposition is shown in Fig.[T5]for the ZUMA runs with no B field (top half), 
and with a uniform Bz =70 T field (bottom half). Table [T] lists the fraction of injected 
hot electron energy deposited in different regions. With no B field, the hot electrons 
propagate essentially freely in the He gas fill. They mostly deposit in the hohlraum wall, 
and a small fraction deposits in the ablator. This is expected based on the solid angle 
subtended by these regions. With a uniform Bz = 70 T, the hot electrons are strongly 
magnetized in the He gas and guided to the capsule. They mostly deposit in the ablator, 
out to a radius comparable to that of the source. This asymmetric preheat, occurring 
mostly in the poles, may drive capsule asymmetries. The energy deposited in the DT 
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Figure 14. Time for 90° angular scatter Tas for: blue curve: 0.96 mgcm“^ of He at Te=l keV 
(typical for hohlraum gas fill in NIF shot N120321), green curve: 1 gcm“^ of CiHi at 200 eV 
(typical for the ablator). Red dashed curves are cyclotron period Tee for Bz — 1 and 70 T. Most 
hot electrons are (strongly, weakly) magnetized in the (He gas, CH ablator). 


layer is ^ 12x higher with the 70 T field. Whether this is a preheat concern depends on 
the spectrum and total energy of hot electrons produced. 

The lack of deposition in high Z with the field means the same hot electron source 
produces many fewer hard x-rays. This is a diagnostics concern, since hard x-rays are 
generally used to deduce hot electrons on NIF. One such principal diagnostic is th e 
FFLEX hard x-ray (> 10 keV) detector ( Dewald et al. 1201 Ol : iHohenberger et al. 20141 ). 
with 10 channels filtered for different energy ranges. Energetic electrons lose energy by 
collisions with background electrons and by bremsstrahlung radiation. Radiation loss 
- EZ X collisional loss, with the two equal in gold for E = 10 MeV. Only electrons 
that deposit energy in high-Z material, such as the hohlraum wall, produce enough hard 
x-rays for EELEX to detect. Hot electrons striking the capsule during the picket have 
been measured on “re-emit” experiments, where the capsule is replaced by a high-Z (e.g. 
bismuth) ball. 


7. Hot electron propagation during peak power: SRSPEAK series 

We now consider ZUMA simulations of propagation of a realistic hot electron source 
produced by Raman scattering on the inner beams during peak power. We use the same 
plasma conditions at 18 ns that were used in the CAPTEST series (i.e., a simulation of 
NIE shot N120321 with the full, incident laser power on each cone), along with conditions 
from a HYDRA run with an initial Bzo = 70 T axial field and the MHD package active. 
The hot electron source has a Maxwell-Jiittner energy spectrum with = 30 keV. 
This temperature is gotten from EELEX data at 18 ns (rise to peak power) on NIE 
shot N130517, which is analogous to N120321 (|Robev et a/.l l2014|) . Once peak power is 
reached, = 18 keV is consistent with EELEX data. The injected angle spectrum is 
dN/dVt = exp[—((6> — 27°)/10°)^], which is directed along roughly the bisector of the two 
NIE inner beams at = 23.5° and 30°. 

We find strong sensitivity to what field lines hot electrons start on - namely, whether or 
not the field lines connect to the capsule. Hot electrons are injected in the three locations 
indicated as sources 1, 2, and 3 in Eig. [5l from r = 0 to 0.06 cm at 2 : = —0.4 cm, and 
from r = 0.12 to 0.18 cm, at 2 : = —0.2 and -0.25 cm. The energy deposition vs. space 
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Figure 15. Hot electron energy deposition for PICKET series of ZUMA runs with no 
magnetic field (top) and a uniform = 70 T field (bottom). 


Deposited energy [J/mnn^ ] per injected hot e- Joule; 1 ns 



-0.4 -0.2 0.0 0.2 0.4 

z [cm] 


Region 

DT gas 

DT layer 
CH ablator 
He gas 

An 

DU 

total 

PICKET, 
no B 

6.56E-5 

2.20E-3 

0.0749 

0.0566 

0.366 

0.428 

0.927 

PICKET, 

B, =70 T 

1.06E-3 (16x) 

0.0261 (12x) 

0.696 (9.3x) 
0.0646 (l.lx) 
4.14E-4 (l.lE-3x) 
4.02E-4 (9.4E-4x) 
0.789 (0.85x) 

SRSPEAK 1, 
no MHD 

4.32E-6 

3.58E-4 

0.406 

0.223 

0.250 

0.0990 

0.979 

SRSPEAK 1, 
B^o = 70 T 

6.26E-5 (14x) 
2.89E-3 (8.1x) 
0.804 (2.0x) 
0.117 (0.52x) 
l.OlE-4 (4.0E-4x) 
1.61E-5 (1.6E-4x) 
0.925 (0.94x) 

Region 

SRSPEAK 2, 

SRSPEAK 2, 

SRSPEAK 3, 

SRSPEAK 3, 


no B 

B^o = 70 T 

no B 

B,o = 70 T 

DT gas 

1.75E-6 

8.95E-9 (5.1E-3x) 

1.44E-6 

5.96E-6 (4.1x) 

DT layer 

1.37E-4 

3.44E-6 (0.025x) 

1.19E-4 

1.26E-3 (llx) 

CH ablator 

0.272 

0.105 (0.39x) 

0.327 

0.576 (1.8x) 

He gas 

0.229 

0.499 (2.2x) 

0.182 

0.248 (1.4x) 

An 

0.335 

0.220 (0.66x) 

0.328 

0.101 (0.31x) 

DU 

0.133 

0.0421 (0.032x) 

0.131 

5.56E-3 (0.042x) 

total 

0.969 

0.866 (0.89x) 

0.968 

0.932 (0.96x) 


Table 1. Fraction of injected hot electron energy deposited in different regions, for PICKET 
(1 ns) and SRSPEAK (18 ns) series of ZUMA runs, (x) is ratio of with H/MHD to without. 
Fractions do not sum to unity because some hot electrons escape from problem boundaries. 


is plotted in Fig. [161 Ihe total into various materials is given in Table [TJ With no 
MHD, the fraction of hot electron energy deposited to DT varies from (1.2 — 3.6) x 10“^ 
over the three sources. The field strongly magnetizes the hot electrons in the He fill gas 
for all three sources, as in the PICKET series of Section |6l Also like the PICKET series 
electrons from source 1, in the LEH, are guided to the capsule. The deposition in DT, CH, 
and He is greatly increased compared to the no-MHD case. For source 2, located deeper 
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in the hohlraum and off axis, electrons are injected on field lines that do not connect to 
the capsule. The resulting deposition in He gas is significantly increased compared to the 
no-MHD case, while that into the CH ablator and especially the DT layer are reduced. 
The situation reverses for source 3, which is slightly closer to the capsule in z than source 
2. Some electrons now start on field lines that connect to the capsule, which results in 
much higher DT deposition. It is not presently known where in the hohlraum SRS hot 
electrons are produced, so we cannot say whether the field increases or decreases DT 
deposition. As with the PICKET series, the fraction of hot electron energy deposited in 
high-Z material is lower with the field, especially for source I. Hard x-ray diagnostics 
may thus not be reliable indicators of hot electron preheat. 


8. Conclusions 

This paper gave results of HYDRA rad-hydro simulations with no MHD effects, and 
with MHD and a 70 T initial axial B field. The field is essentially frozen-in to the highly 
conductive plasma, and gets advected with the radial expansion of the capsule and wall. 
This results in field lines that roughly follow contours of ablated material. The magnetic 
pressure is much less than material pressure. The principal hydro effect of the field is 
reduced electron heat conduction perpendicular to it. This gives a hotter hohlraum fill, 
especially in gold, and a wider channel between the capsule and equator wall. Less inner- 
beam absorption occurs before they reach the wall, which increases the equatorial x-ray 
drive. Inner-beam Raman scattering may be reduced by the hotter fill, in addition to the 
lower power needed to achieve a round implosion. 

We also presented hot electron propagation studies with ZUMA, using plasma condi¬ 
tions from HYDRA. Mono-energetic test cases with plasma conditions from early peak 
power (18 ns) show a minimum hot electron energy of 125 keV incident on the capsule 
is required to reach the DT layer. The energy coupled to the layer maximizes at 13% for 
185 keV electrons, and drops with energy above that. Using plasma conditions during the 
early-time picket (I ns) with no field, we find a small fraction (2 x 10“^) of hot electron 
energy from a two-plasmon-decay relevant source couples to DT. With a uniform 70 T 
axial B field, the hot electrons are magnetized in the He fill gas, guided to the capsule, 
and the DT coupling increases by I2x. This may not be a preheat concern, since picket 
pulse shaping has been shown on NIF to significantly decrease the hot electron source. 
ZUMA simulations using plasma conditions at 18 ns, with a source motivated by inner- 
beam SRS, show an imposed field can greatly increase or decrease hot-electron coupling 
to DT. This depends on whether electrons are produced on field lines that connect to 
the capsule. 

Imposed magnetic fields may enhance hohlraum performance by improving inner-beam 
propagation and reducing Raman scattering during peak laser power. This is in addition 
to the primary benefit of reducing electron-heat and alpha-particle loss from the hotspot. 
One concern is possible increase in DT fuel preheat due to the field guiding hot electrons 
to the capsule. Work is underway on a pulsed-power field generator for NIF, and we look 
forward to hohlraum experiments in the next few years. 

We gratefully acknowledge fruitful conversations with H. F. Robey, J. D. Salmonson, C. 
A. Thomas, J. Hammer, and D. E. Hinkel. This work was performed under the auspices 
of the U.S. Department of Energy by Lawrence Livermore National Laboratory under 
Contract DE-AC52-07NA27344. Partly supported by LLNL LDRD project I4-ERD-028. 
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Figure 16. Hot electron energy deposition for SRSPEAK series of ZUMA runs without MHD 
(top r > 0) and with MHD (bottom r < 0). Top, middle, and bottoms plots are for source 
locations 1, 2, 3 indicated in Fig. [5l 
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